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Abstract 



Examples of the construction of Hamiltonian structures for dynamical sys- 

o ' 

terns in field theory (including one reputedly non-Hamiltonian problem) with- 



out using Lagrangians, are presented. The recently developed method used 
requires the knowledge of one constant of the motion of the system under 
consideration and one solution of the symmetry equation. 

I. INTRODUCTION 



Hamiltonian theories have been widely used in almost all areas of physics. The usual 
approach consists in constructing the momenta (canonically conjugate to the coordinates) 
and the Hamiltonian starting from a Lagrangian formulation of the system under consider- 
ation. Nevertheless, there has been increasing interest in studying non standard procedures 
to produce Hamiltonian structures starting from the equations of motion only, without us- 
ing a Lagrangian ||. These approaches usually deal with systems which are naturally 
described in terms of non-canonical variables and the Hamiltonian formulation is conse- 
quently written using non-canonical Poisson matrices, which in many cases are singular. 
Due to the lack of a general procedure to generate Hamiltonian structures from scratch, in 
most instances the results have been obtained by extremely inspired guesswork. 

The purpose of this article is to present examples of the application of a newly devised 



method [|10l , |TT]l which constitutes a general technique for the construction of Hamiltonian 



structures for dynamical systems. We discuss systems described by non-linear equations in 
field theory and also some linear equations such as the time dependent Schrodinger equations 
and the heat equation, which according to the usual belief, is clearly non-Hamiltonian ||. 

In Section 2 we give a brief outline of the method and the rest of the Sections are devoted 
to the construction of the examples. 



II. BRIEF OUTLINE OF THE METHOD 



The contents of this Section summarize partially the results of |IU |. We use classical 
mechanical notation nevertheless, the results may be easily applied to field theory as well, 
as we will see at the end of this section. 

Consider a dynamical system defined by equations which have been cast in first order 
form, 

dx a 



/V) a, 6 = 1, ,N. (1) 



dt 

A Hamiltonian structure for it consist of an antisymmetric matrix, J ab (x c ) and a Hamil- 
tonian H(x c ) such that J ab is the Poisson bracket for the variables x a and x b (which are 
non-canonical in general) and H is the Hamiltonian for system ([!]). In addition to its an- 
tisymmetry, the matrix J ab is required to satisfy the Jacobi identity and to reproduce, in 
conjunction with the Hamiltonian H, the dynamical equations (p, i.e., 

r\ d j dc + j bc , d j da + r a , d J db = o , (2) 

and, 

^ = r- (3) 

It has been proved |10] that one solution to the problem of finding a Hamiltonian structure 
for a given dynamical system is provided by one constant of the motion which may be used 
as the Hamiltonian H, and a symmetry vector r] a which allows for the construction of a 
Poisson matrix J ab . The constant of the motion and the symmetry vector satisfy, 



C f H = 0, (4) 
(d t + C f )rf = , (5) 



respectively, where Cf is the Lie derivative along / (for a definition, see [12|], for instance). 
In addition, it is required that the deformation K of H along r) a , 

K = ^ a = ' (6) 

be non-vanishing. The Poisson matrix J ab is constructed as the antisymmetrized product of 
the flow vector f a and the "normalized" symmetry vector rf/K, 

J ab = ^(f a V b ~ fV) ■ (7) 

The Poisson matrix so constructed has rank 2 and it is, therefore, singular in many 
instances. Adding together two Poisson matrices constructed according to (|7|) will not 
increase its rank. It will just redefine the symmetry vector used to construct it. One 



method to increase the rank of such a Poisson matrix is presented in |T0 . 

Let us now deal with systems with infinitely many degrees of freedom defined by some 
field <f)(x,t), where x denotes the coordinates of a point in space. The dynamical equation 
(0) now is, 

j>{x,t) = F[<f>,x) , (8) 

where F is a functional of (f) for every point x in space. All the above discussion remains 
valid replacing x a by 4>(x), tensors T ab '" c (x a ) by functionals which depend on some spatial 
coordinates, Q[<f>, x, y . . . , z] and partial derivatives d/dx a by functional derivatives 5/5<j)(x). 
Details are given in [[|. 

III. HEAT EQUATION 

Let us consider the heat equation, 

u t = u xx . (9) 
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It is easy to see that 

r)(x) = A , (10) 

is a symmetry transformation for it, where e is any real number. 

If we assume periodic boundary conditions in x G [—a/2, a/2] then the quantity 

pa/2 

H= udx (11) 

J -a/2 

is conserved. 

The deformation of H along 77 is given by 

,a/2 $H pa/2 

C V H = / dx 1 — r ^T](x) = A dx = Aa , (12) 

J-a/2 du{X) J-a/2 

where all the integrals are taken from —a/2 to a/2. 

We can now construct a Poisson matrix associated with the Hamiltonian ([11]) for equation 



J(x,y) = Uxx ^ UM . (13) 

(X 

Let us explicitly show that J ab together with H provide a Hamiltonian formulation for 
the heat equation. In fact, 

H 2 5H 
[u(x),H] = / dyJ(x,y)- 



a/2 ' Su(y) 

1 ra/2 I ra/2 

= -u xx (x) / dy / u yy (y)dy 

a J-a/2 a J-a/2 ' 

So we see that, 

ii=[u,H], (14) 

as required. It is worth noting that according to folk tradition, this equation cannot be 
endowed with a Hamiltonian structure. In Salmon's words ||: "By anyone's definition, 
(0.1) is non-Hamiltonian" . In his paper, (0.1) is the heat equation subject to periodic 
boundary conditions. 



IV. TIME DEPENDENT SCHRODINGER EQUATIONS 



Consider the following equations of motion 

il>t = iWx* + V(x,t)il>) = f, (15) 

and its complex conjugate 

1>% = -iW„ + V(x,tW)=f . (16) 

Note that we discuss the case of a time dependent potential. We use one dimensional 
notation for simplicity only, our results hold irrespective of the dimensionality of space. 

It is a straightforward matter to realize that multiplication of the variables ip and ip* 
by two different constants, (1 + A) and (1 + fi\) respectively, constitutes a symmetry trans- 
formation for the Schrodinger equations ( |T5| ) and (pp. The infinitesimal version of this 
transformation is 

n = i> , (17) 

and 

rf = pi/;* . (18) 
The usual probability conservation statement means that 

H = J ip*(x)ip{x)dx , (19) 
is a conserved quantity as it may be easily proved. Its deformation K along the symmetry 



vector defined by fll7[) and fll8|) may be written in terms of variational derivatives as 



K ^ I m i,{x)dx + f ^' nx)dx • (20) 

and a straightforward calculation yields 

K=(l + n)H^0, (21) 



for any /x such that 1 + // ^ . The Poisson structure is defined by 

hK*),V(v)] = -^(/(x)77(y)-77(x)/(y)) , (22) 

[^(x),^(y)] = ^(/(i)ij*d/) - /*(*)) , (23) 

and 

ir(x),r(y)] = ^(r(oc) V *(y)- V *(x)f*(y)) . (24) 

We have thus constructed a family of Poisson structures which depend on the parameter 
fi for a given Hamiltonian H. Note that in spite of the time dependent character of the 
Schrodinger equations (|i~5|) and ( |i~6l) the Hamiltonian ( |i~9|) is conserved, while the usual 
Hamiltonian contains the time dependent potential V(x, t) and it is not conserved. On the 
other hand, our Poisson structures are time dependent. This is, as far as we know, a novel 
feature for Poisson matrices, which means that, even in the case of a regular matrix, the 
Hamiltonian structure is not derivable from a Lagrangian. 



V. KORTEWEG DE VRIES EQUATION 

The equation of motion is 

u t = -uu x - u xxx = f . (25) 
It is not difficult to see that symmetry transformation for it is given by [|Kj 



77 = (-2m - xu x + 3t(uu x + u xxx )) . (26) 

In fact, to prove that ( p6|) is a symmetry transformation for the KdV equation, it is enough 
to check that 77 satisfies the symmetry equation (|5]) with / defined by (f2~5|). 

To get a Hamiltonian structure for the KdV equation we need to construct constants of 



the motion which be non trivially deformed by 77. In [10] the energy 
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Hi= J u 2 dx , (27) 

was used as a constant of motion to complete the Hamiltonian structure. Note that the 
deformation K\ of H\ along rj is non-vanishing 



/5H 
—^ri^dx = -3#i . (28) 
du[x) 



The Poisson structure Ji(x,y) is given by 



Ji(x, y) = ^(f(x)rj(y) - f(y) V (x)) . (29) 



Consider now Ho 



u 2 

H 2 = l(-^-+u 3 )dx, (30) 



which is also conserved. Its deformation K 2 along rj is 



K 2 = / — -n(x)dx = -5H 2 . (31) 
J ou(x) 

The Poisson structure J 2 (x,y) is given by 

J 2 (x, y) = -L(f( x )r,(y) - f(y) V (x) ) . (32) 

We have been able to construct two different Hamiltonian structures based on one symme- 
try vector rj and two conserved quantities Hi and H 2 as the Hamiltonians of each of the 
structures. The Poisson structures J\ and J 2 are built as the antisymmetric product of the 
evolution vector / and the symmetry vector rj normalized using the deformations K\ and K 2 
of H\ and H 2 , respectively. Note that a different but closely related scheme exists in which 
no normalization of the antisymmetric product is needed provided that the Hamiltonians 
— (1/3) log if i and —(1/5) log if 2 be used instead of Hi and H 2 respectively. In this case the 
Poisson matrix for both Hamiltonian structures is exactly the same one. 

A similar construction may be performed with the rest of the constants of the motion 
which belong to this family, as they appear, for instance, in |0. We have then constructed 



a set of infinitely many Hamiltonian structures for the KdV based on one symmetry trans- 
formation and different constants of the motion. Let us remark that as it was mentioned in 
Section II, the matrices J(x, y) constructed according to (0) have rank 2 hence, there are 
many Casimir functions which have vanishing Poisson bracket relations with any other dy- 
namical quantity. A few words regarding the construction of Casimir functions for Poisson 
structures such as these seem in order. Time independent Casimir functions are, of course, 
constants of the motion. Consider now a different evolution along another parameter (call 
it s) which is given in terms of the symmetry vector r] a . In other words, consider 

dx a 



ds 



rj a (x b ) a, 6=1, ,N . (33) 



Deformation of constants of the motion along rj a may be viewed as the evolution of such 
constants in the parameter s. A Casimir function is such that its evolutions (both in time 
and in the s parameter) vanish. So, the construction of Casimir functions may be viewed as 
the search of entities which are simultaneously constant for both the time and s evolutions. 
To illustrate this, take for instance Eqs. (p8|) and (|3l| ) and rewrite them as 

dHt 



ds 

and 

dH, 



-3#i . (34) 



-5H 2 . (35) 



ds 

solve them and eliminate the parameter s to get that Hi 1 ' 3 H 2 ~ 1 ^ 5 is a Casimir function for 
both of the Poisson matrices ( f29|) and ([32|). Similarly, one can get Casimir functions for 
other Poisson matrices. 

VI. BURGERS EQUATION 

Consider the Burgers equation, 

u t = u xx + 2uu x = f , (36) 



S 



x being in [—a/2, a/2]. It is straightforward to prove that 

r a/2 



Hi 
H 2 



-a/2 
a/2 

-a/2 



dxu(x) 



dxexp 



a/2 



u{y)dy 



(37) 
(38) 



are conserved quantities for it, if we assume that the field vanishes at the boundary ±a/2. 
A symmetry transformation for Eq. ([36]) is given by 



rj(x) = u(x) exp 
To see this, it is enough to show that 



u 



a/2 



In fact, 



fa/2 f 6 f(x) 5 V (X) " 



dy^Y^Viy) = / dy——-f{y) 

-a/2 0U{y) J -a/2 0U{yJ 



-a/2 

The deformation of Hi along rj does not vanish, 

r«/2 , SHi 



exp ( J u(y)dy j (u xx + uu x - w 3 ) . 



raj 

Ki = C v Hi = / 

J —a 



dx- 



T](x) 



/2 5u(x) 

a/2 

dxu(x) exp 

-a/2 

1 _ exp(-Hi) . 



a/2 



dwu(w) 



It may be proved that the deformation of H 2 along r\ is also non vanishing, 
K 2 = C V H 2 -- 



a / 2 , SH 2 
ax- — ^-r](x 



-a/2 5u(x) 
a/2 / rz \ 

dzexp / u(z) 

-a/2 \J-a/2 J 

Ho- a. 



a/2 



do; exp 



a/2 



dwu(w) 



(39) 



Therefore, we may construct Hamiltonian theories for the Burgers equation using either Hi 
or H 2 as Hamiltonians. The appropriate Poisson matrices are 



and 



Mx,y) = f{x)v{y) - f{y)ri{x) , (40) 



J2{x,y) = , (41) 

&2 



respectively. 

It is easy to check, for example, that, 



e Hl - 1 

C = 7T^' (42) 

n 2 — a 



is a Casimir for both of the Poisson brackets defined by (fiTf ) and (41), where we have used 
the approach described at the end of the preceding Section. 

VII. HARRY DYM EQUATION 

The Harry-Dym equation is 

u t = in- 1 ' 2 ) = f , (43) 

\ / XXX 

where x G [—a/2, a/2]. If we assume periodic boundary conditions, H\ and H 2 

/•a/2 

Hi— udx , (44) 

J — a/2 
pa/2 

H 2 = / u 1/2 dx , (45) 

J -a/2 

are conserved quantities. Note that a may be set equal to oo. 
Let us define the vector field 

£(x) = Au — Bxu x , (46) 
where A and B are real constants. Let us now compute the Lie derivative of it along /, 

£,£ = 3 (i> + *)(«-•/»)_ = 3 (lA + *)/. (47) 

Therefore, rj 
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T, = -3t(±A + B)f + £, (48) 

is a symmetry transformation and can be used to construct a Hamiltonian theory provided 
it deforms non trivially some Hamiltonian. 

This is exactly the case for Hi and H 2 . In fact, 

Ki = C v Hi = (A + B)Hi , (49) 
K 2 = C V H 2 =(^ + B^H 2 , (50) 

so we have one family of Poisson matrices associated with Hi and another one associated 
with H 2 . They are 

( vT 1 / 2 ) (Au(y) - Byu v ) - ( vr 1 ' 2 ) (Au(x) - Bxu x ) 
Ji(x, y) = ± ^ Jjm , (51) 

and 

(V 1 / 2 ) (Au(y) - Byu v ) - (u- 1 / 2 ) (Au(x) - Bxu x ) 
J 2 (x, y) = ± ^ -A ^ , (52) 

respectively. Of course, we must be careful to choose A and B in such a way that either Ki 
or K 2 (or both) be non-vanishing. 



VIII. CONCLUSIONS 

We have constructed Hamiltonian structures, without using Lagrangians, for several 
linear and non-linear systems of partial differential equations (field theory) based on a 



recently devised method [|iO[|iT| , which needs only of the knowledge of a constant of motion 
and a solution of the symmetry equation. The examples include the heat equation which has, 
up to now, been considered to be non-Hamiltonian |J. The structures found are singular 
in the sense that there exist Casimir functions (which have vanishing Poisson bracket with 
any dynamical variable). This feature is present in many other Hamiltonian theories, and 
it is sometimes unavoidable (as it is the case of gauge and constrained systems), and it 
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also appears in Hamiltonian descriptions of some fluids. We should stress that being able 
to produce a Hamiltonian structure (albeit singular) for a system of differential equations 
constitutes progress with respect to the situation of having no Hamiltonian structure at all. 
More examples will be discussed in forthcoming articles. 
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